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Abstract 

A new and powerful mean field scheme is presented. It maps to 
a one-dimensional finite closed chain in an external field. The chain 
size accounts for lattice topologies. Moreover lattice connectivity is 
rescaled according to the GM law recently obtained in percolation 
theory. The associated self-consistent mean-field equation of state 
yields critical temperatures which are within a few percent of exact 
estimates. Results are obtained for a large variety of lattices and di- 
mensions. The Ising lower critical dimension for the onset of phase 
transitions is = 1 + |- For the Ising hypercube it becomes the 

Golden number di = ^'^^ ■ The scheme recovers the exact result of 
no long range order for non-zero temperature Ising triangular antifer- 
romagnets. 



* Laboratoire associe au CNRS (UMR 7603) 



1 Introduction 



For many years basic mean-field theory has been apphed to a huge variety of 

problems. It is a very simple way to tackle collective phenomena [1]. Most of 
the time it yieds a correct qualitative description. However, quantitatively 
the results are very poor. In particular, all aspects of the critical behavior 
are .grossly misrepresented [2]. Moreover some results are even wrong like 
for instance the existence of long range order for the Ising system at both 
one dimension for ferromagnets and at two dimensions for the triangular 
ant if erromagnet s . 

The Bethe scheme [3] was then introduced to extend the Weiss one spin 
approach [4] to a cluster of fiuctuating spins. For decades it has been looked 
upon as a solid imrovement. It yields no long range order at one dimension 
and even becomes exact for the hypercube at infinite dimensions [1] . However 
quantitative results are yet rather poor. Moreover it was demonstrated re- 
cently that the Bethe scheme violates systematically translational invariance 
[5]. As such it is forbidden by symmetry. 

In this paper we present a new and powerful mean field scheme. It em- 
bodies the Bethe idea of including a fiuctuating spin cluster yet preserving 
the overall lattice translational invariance. In addition, the connectivity be- 
tween fluctuating clusters is rescaled according to the GM law introduced 
few years ago in percolation theory [6]. 

Associated critical temperatures are calculated for a large variety of lat- 
tices and at several dimensions. Discrepancies with available exact estimates 
are only within few percent. The lower critical dimension for the onset of 
phase transitions is found to be = 1 + ^ for Ising systems. It turns to 

he Golden number di = for hypercubes (q = 2d). In the case of the 

triangular Ising antiferromagnet the exact result of no long range order is 
reproduced [7] 

2 Revisiting the Bethe aprroximation 

Mean-field theory is a one-site approach which first breaks the lattice symme- 
try by discriminating bewteen fluctuating degrees of freedom and averaged 
ones [1]. Two interpenetrated lattices are thus deflned. Equating the thermal 
average of fiuctuating degree of freedom to the already averaged ones restores 
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the inital lattice symmetry. Simultaneously a self-consistent equation of state 
is obtained. 

To implement a Bethe scheme [3] on a lattice 3 distinct interpenetrated 
sublattices (A, B, C) must be introduced. First the fluctuating center (A), 
then the fluctuating nearest neigbhors (B) and last the mean field nearest 
ncigbhors (nn) of the nn (B) not including the center (A). From the A-spin 
plus these B and C shells, a cell is constituted to pave the whole space and 
reproduce the full lattice topology. 

Cluster center (A) has thus all its nn spins (B) which are fluctuating while 
surface cluster spins (B) have mean-field nn spins (C) and one nn fluctuating 
spin (A). Simultaneously mean-field spins (C) have all their nn which are 
fluctuating spins, making their environnement identical to the cluster center. 

On this basis the Bethe requiremnt < Sa >—< Sb > is not compatible 
with the equality mc =< Sb > which should also hold to ensure translational 
invariance. The Bethe topology is therefore forbidden by symmetry. It is not 
the case for one-site Weiss theory. For a detailled demonstration see [5] 

Last but not least, it is worth noticing it is indeed this very symmetry 
problem which makes the Bethe approach exact on the Caley tree lattice. 
This lattice does not exhibit translational invariance by construction. This 
symmetry breaking was overlooked for several decades. 

3 A new powerful mean field schme 

From the discovery of a systematic Bethe induced symmetry breaking arises 
the question of the possibility to indeed extend a mean field treatment to 
more than one site. 

Above analysis of the Bethe scheme emphazises the role of the cluster 
center in the irreversible breaking of the symmetry. It hints to avoid such 
a fluctuating center. One way to achieve this constraint is to use compact 
closed linear loops within the lattice topology. For instance compact 4-spin 
squares and 3-spin triangles for respectively square and triangular lattices. 

Each one of these plaquettes is then set respectively as A-species (fluc- 
tuating) and B-species (mean field) with a staggercd-like coverage pattern. 
A-plaquettes (B-plaquettes) have thus all their nn plaquettes as B-plaquettes 
(A-plaquettes). 

For a given plaquette, each spin has two nn spins of the same species 
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within the plaquette itself and {q — 2) nn spins of the other species belonging 
to nn plaquettes. At this stage we have a series of fluctuating one-dimensional 
closed chains in an external field h. The number of spins N in each chain is 
determined from lattice topology. It is = 4, = 3, = 3 and A'" = 6 for 
respectively square, triangular, Kagome and Honeycomb lattices. 

it is the interactions with nn mean-field spin plaquettes which produce 
the field h. We have h = 6Jm where 6 accounts for connectivety to B- 
sublattices, J is the nn coupling constant and m the averaged magnetization 
on the B-sublattice. The problem can now be solved exactly. In particular, 
the chain site magnetization is [1], 

<5.>^^exp2/f{ |;;;^w:) >.. (1) 

at order one in h. Here i G A-plaquettes. K = where /c^ is the Boltzman 
constant and T the temperature. 

Putting < Si >— m restores the initial lattice symmetry. It is indeed 
possible since only two sublattices were involved which was not the case for 
the 3 sublattice Bethe scheme. The self-consistent equation of state is, 

^^.X.xp2A-{ |;;;^ff) }^ + .... (2) 

at order one in m and using h — SJm. To solve Eq. (2) needs to determine 
the value of S. 

It is then worth to evoke a recent work on percolation thresholds, the 
GM law [6]. It shows that relevant connectivity variables for site and bond 
dilution are repectively {d — l){q— 1) and i^tlll^iili in other words, for site 
percolation, the number of possible directions (g — 1) from a given site, has 
to be multiplied by {d — 1). For bond percolation this effective number of 
site directions has to be divided by dimension d. Using these variables, the 
GM law was found to yield all percolation thresholds for all Bravais lattices 
at all dimensions [6]. 

This percolation finding suggests to consider here a rescaled connectivity 
between closed loops instead oiS — q — 2. Using above counting, we first start 
with q instead of (g— 1) since now dealing with pair exchange interactions and 
not percolation. Second we renormalize q by (d—l) giving q{d— 1). However 
the 2 neighboring sites which are treated exactly within the closed loop have 
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to be substracted from the effective number of sites which gives q{d— 1) — 2. 
Moreover, interactions being related to bonds, we divide this number by d 
as for bond percolation. These considerations lead to a connectivity. 



Results 



We can now check the validity of our simple symmetry preserving model with 
respect to critical temperatures. From Eq. (2) we get, 

5i^«exp2K«|il^^^^^i^^2^|-l (4) 
OA, exp I ^ tanh(i^G)iV) > - ^ • l^j 

The trivial connectivity counting 5 — q — 2 already improves Weiss model. 
For instance = 0.29 in the sqiiare case and Tc — a,t d — 1. We now 
proceed using Eq. (3) for connect ivty. 

For the square case (g = 4, N = 4), 5 = 1 which gives = 0.4399. 
Exact result is — 0.4407. In the case of triangular lattice (g = 6, N = 3), 
= 0.2919 with 5 = 2 while the exact estimate is = 0.2746. For 
Kagome (g = 4, N = 3) S = 1 yielding — 0.4649 for an exact estimate of 

= 0.4666. And = 0.6160 for the honeycomb lattice (g = 3, = 6) 
where 5 = | for an exact estimate of = 0.6585 (see Table I). 

Going to d = 3 imposes to restrict the plaquette size to A" = 4 since a 
one-dimensional loop cannot embody a three-dimensional topology. However 
there exits a (i— dependence through Eq. (3). We get 5 = y, 5 = y and 
5 = Y for respectively cubic, fee and bcc lattices. Corresponding critical 
temperatures are given by = 0.2012, 0.1568, 0.1096 respectively for 
exact estimates of 0.2217, 0.1575, and 0.1021 (see Table I). 

Critical temperature estimates [8, 9] are available for the hypercube at 
d = 5, 6, 7. These are = 0.1139, 0.0923, 0.0777 respectively. 

To get the d ^ oo asymptotic limit of our model we take both q ^ oo 
and J under the constraint qj = est. From Eq. (3) connectivity limit 
is 5 — > g(l — ^) which gives always, 

5^q (5) 
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at leading order. Indeed q diverges always quicker than ^.even for /cc- lattices 
where q = 2d{d — 1). In turn Eq. (4) becomes, 



K^-^, (6) 

which is the mean- field result [1] as expected in the d ^ oo limit. 

To evaluate the sensibility on the loop size, it is fruitful to expand Eq. 
(4) in powers of K. It gives 

(1 + 2K^ + ... + - {K^r + -Oil - {K^r + •••) = J • (7) 

Since A/" > 3, a simple analytic expression is obtained only at order one, 

K^-\- (8) 

At two dimensions Eq. (8) gives K'^ = 1, |, 1, 2 for respectively the square, 
triangular, Kagome and for honeycomb lattices. These results are rather 
poor and shows the importance of the finite value of N which embodies part 
of the lattice topology. 



5 Conclusion 

Wc have presented a very simple self-consistent model which yields rather 
good values for critical temperatures within a few percent of exact results. 
Besides a rescaled lattice connectivity, the finite length of the loops is also 
taken into account. This new scheme represents a substantial improvement 
over existing mean-field cluster approximations. 

We can also determine from our model a lower critical dimension for 
phase transitions. It comes from the condition /i = for which we have a 
one-dimensional finite system. Such a system has no long range order at 
T ^ 0. Phase transitions are thus obtained only in the range h ^ which 
gives q{d — l)>2 leading to, 

di = l + - . (9) 
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For the Ising hypcrcubc {q = 2d) it becomes the Golden number di = ' 
which excludes the d = 1 case and contains d = 2 as it should be. Last but 
not least, applying our scheme to the Ising triangular antiferromagnet[7] we 
do revover the exact result of no long range order at non-zero temperatures; 
contrary to usual mean field approaches. 
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Dimension 


Lattiee 


Q 


5 






= 2 


Square 


4 


1 


0.4407 


0.4399 




Honeyeomb 


3 


\ 

2 


0.6585 


0.6160 




Triangular 


6 


2 


0.2746 


0.2837 


d = 2 


Kagome* 


4 


1 


0.4666 


0.4649 


d = 3 


Diamond 


4 


2 


0.3698 


0.2857 




se 


6 


10 

3 


0.2216 


0.2012 




bee 


8 


14 

3 


0.1574 


0.1568 




fee 


12 


22 
3 


0.1021 


0.1096 


d^A 


se 


8 


22 
4 


0.1497 


0.1380 




fee 


24 


23 
2 




0.0749 


d^h 


se 


10 


38 
5 


0.1139 


0.1064 




fee 


40 


158 

5 




0.0298 


= 6 


se 


12 


29 
3 


0.0923 


0.0869 


d^7 


se 


14 


82 
7 


0.0777 


0.0737 



Table 1: from this work compared to "exact estimates" Kl taken from [8, 
9]- 
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